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Quasi-One-Dimensional Modeling of 
Pulse Detonation Rocket Engines 


C. I. Morris* 

NASA/George C. Marshall Space Flight Center 
TD40/Propulsion Research Center 
Marshall Space Flight Center, AL USA 35812 

Pulse detonation rocket engines (PDREs) offer potential performance improvements over conventional de- 
signs, but represent a challenging modeling task. A quasi l-D, finite-rate chemistry CFD model for a PDRE is 
described and implemented. A parametric study of the effect of blowdown pressure ratio on the performance 
of an optimized, fixed PDR£ nozzle configuration is reported. The results are compared to a steady-state rocket 
system using similar modeling assumptions. 


Introduction 

P ULSED detonation rocket engines (PDREs) have gen- 
erated considerable research interest in recent years ^ 
as a chemical propulsion system potentially offering im- 
proved performance and reduced complexity compared 
to conventional rocket engines. The detonative mode of 
combustion employed by these devices offers a thermo- 
dynamic advantage over the constant-pressure deflagra- 
tive combustion mode used in conventional rocket engines 
and gas turbines. However, while this theoretical advan- 
tage has spurred a great deal of interest in building PDRE 
devices, the unsteady blowdown process intrinsic to the 
PDRE has made realistic estimates of the actual propul- 
sive performance problematic. The recent review article by 
Kailasanath^ highlights some of the difficulties in compar- 
ing the available experimental measurements with numeri- 
cal models. 

In a previous paper by the author,^ parametric studies of 
the performance of a single, straight-tube PDRE were re- 
ported. A l-D, unsteady method of characteristics code, 
employing a constant-y assumption behind the detonation 
front, was developed for that study. Models of this type 
are computationally inexpensive, and are particularly use- j 
ful for parametric performance comparisons. For example, 
a plot showing the specific impulse of various PDRE and 
steady-state rocket engine (SSRE) configurations as a func- 
tion of blowdown pressure ratio (Pinitiai/Pfinai) is shown 
in Figure 1. Note that the SSRE performance calcula- 
tions employ the same constant- y assumption used in the 
PDRE calculations. The performance curves indicate that 
a straight-tube PDRE provides superior specific impulse, 
compared to a SSRE with a sonic nozzle, over the entire 
range of pressure ratios. Note, however, that a straight- 
tube PDRE in general does not compare favorably to a 
SSRE fitted with an optimized converging-diverging (C-D) 
supersonic nozzle, particularly at the high pressure ratios 
typical for boost or in-space rocket applications. This result 
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is largely due to the choked outflow from a straight-tube 
PDRE. However, calculations of an ideal expansion of the 
PDRE outflow show that if a dynamically optimized super- 
sonic nozzle could be fitted to a PDRE, then the specific 
impulse of the device would exceed that of a comparable 
SSRE. While such a nozzle is a considerable idealization, 
it is clear that nozzle design and optimization will play a 
critical role in whether the performance potential of PDREs 
can be effectively realized in practice. 

The purpose of this paper is to report efforts to study 
the effect of nozzles on PDRE gasdynamics and perfor- 
mance. Details of the quasi l-D, finite-rate chemistry CFD 
model developed by the author are provided first. A para- 
metric study of the effect of blowdown pressure ratio on the 
performance of an optimized, fixed PDRE nozzle configu- 
ration is then reported. The results are then compared to a 
SSRE system using similar modeling assumptions. 



Fig. 1 Performance comparison of a pulse detonation rocket 
engine (PDRE) with a conventional steady-state rocket engine 
(SSRE) equipped with both sonic and optimized supersonic 
converging-diverging nozzles. The specific heat ratio, y, is 
held constant in all models. The final blowdown pressure at 
the end-wall in the PDRE is equal to the ambient pressure 
(Pfinal * Fambient)* Propellant mixture: stoichiometric H2-02. 
Propellant initial conditions: Pinitial = ^ Rl™* Tiniual = 300 K. 
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Governing Equations 

The quasi-one-dimensional Euler equations are the gov- 
erning model used in this study. A quasi- 1-D model is at 
best an approximation of a pulsed detonation or steady- 
state rocket engine. Real detonation waves exhibit clear 
multi-dimensional behavior. Additionally, flow separation 
in rocket engine nozzles is strongly dependent on multi- 
dimensional and viscous effects . However, as the primary 
goal of this study is to determine the effect of nozzles on 
pulse detonation rocket systems to first order, the quasi- 1- 
D Euler equations were deemed a reasonable model. The 
time-dependent form of these equations, including the ef- 
fect of finite-rate chemistry, is written as follows: 


W 1 3FS 
dt S~^ 


H + W 


(I) 


where 5 represents the cross sectional area of the duct, 
and is a function of distance, x. The state vector, U, and 
convective flux vector, F, are given by 
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The source vectors accounting for the effect of area change, 
H, and for the effect of finite-rate chemistry, W, are 
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Here, pi , . . . , Pns represent the mass densities of the ns 
individual chemical species in the problem. The total gas 
density p is sum of the individual species densities, p = 
5Z!t=i ^ velocity in the x direction, p is the 

gas pressure determined by the ideal gas law. 


p = y^RT 


(4) 


k^\ 


where Mk is the molecular weight of the species, R is the 
universal gas constant, and T is the gas temperature, e is 
the total energy per unit mass of the gas, ^ € -h 

where € is the specific energy per unit mass, and for the 
thermally perfect gas model considered here f is a function 
of temperature and gas composition, 




(5) 


k=\ 


where Yk represents the mass fraction of species k, = 
Pk/py and €k{T) is the specific energy per unit mass of 


species k, as a function of temperature. In this work, the 
ideal gas thermodynamic fits of McBride et al."^ are used to 
calculate the energy, specific heat, and Gibbs free energy 
for the 9 species in the problem (N 2 , O 2 , H 2 , OH, H 2 O, H, 
O, HO 2 , and H 2 O 2 are included). 

Chemistry Model 

The source terms for finite-rate chemistry, Wk in Eq. 3, 
are calculated from the sum of relevant chemical reactions 
for each species. Thus, for a chemical mechanism of nr 
elementary reactions, each reaction is expressed as 




( 6 ) 
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where the v[ ^ are the stoichiometric coefficients of species 
/ on the reactant side of reaction /i, while the are the 
corresponding stoichiometric coefficients on the product 
side. The production term for each species may be writ- 
ten as 

^Pk 


m^ — = Mk - <n) 


n=l 


(7) 


/=! /=! 
and k / and kh represent the forward and backward reaction 
rates, respectively. Details of the chemical kinetics mech- 
anism for H 2 -O 2 combustion used in this work are shown 
in Table 1. This mechanism was developed by Petersen 
and Hanson^ to model H 2 -O 2 ignition at the elevated pres- 
sures typical of practical high-speed propulsion systems. 
The rate coefficients for reactions without pressure depen- 
dence take on the conventional Arrhenius form: 


kf{T) = AT^ exp(-Eact//?7’) 


( 8 ) 


where A is a constant, and Eact is the activation energy in 
cal/mol. All reactions in the mechanism are reversible. The 
reverse reaction rate, kh, is calculated form the forward rate 
(Eq. 8) and the equilibrium constant. 

Pressure-dependent reactions are modeled using the ap- 
proach given by Kee:^ 


kf=koo[Pr/{\^Pr)]F 
where the reduced pressure, Pr, is 
^ __ ko[M] 

^00 

and the correction factor, F, is in the Troe^ form: 

1-1 


InF 


1 + 


\nPr -\-c 


-1 

In i 
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( 10 ) 


( 11 ) 


n — d(\n Pr-\-c) 

The Troe centering parameter, Fc is given by 
F, = (1 -a)exp(-F/r**) + aexp(-7/^*)^- 

exp(-rV7') (12) 

where the constants c = —0.4 — 0.67 In Fr, n = —0.75 — 
1.27 In Fo and = 0.14. 
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Table 1 Chemical kinetics mechanism for H2-O2 combustion, from Petersen and Hanson.^ 


Number 

Reaction” 


aI> 

n» 

^act 

1 

0 -|-H 2 ^H + 0 H 


5.00 X 10* 

2.70 

6290 

2^ 

H + O2 -f M ^ HO2 M 


2.80 X 10'* 

-0.90 

0 

3 

H - 1 - O2 + O2 ^ HO2 - 1 - O2 


3.00 X 10“ 

-1.70 

0 

4 

H - 1 - O2 + H2O HO2 - 1 - H2O 


9.38 X 10'* 

-0.80 

0 

5 

H 4 - O2 + N2 ^ HO2 -H N2 


2.60 X 10'* 

-1.20 

0 

6 

H + O2 ^ 0 4 - OH 


8.30 X 10'^ 

0.00 

14413 

7 

H~\-H 02^02 + H 2 


2.80 X I0'3 

0.00 

1068 

8 

H + HO2 ^ OH + OH 


1.34 X I0''‘ 

0.00 

635 

9 

H 4 - H2O2 ^ HO2 4 - H2 


1.21 X 10’ 

2.00 

5200 

10 

0 H-FH 2 :?^H 20 -FH 


2.16 X 10* 

1.50 

5200 

ll^.e 

H2O2 M 

^00 

7.40 X 10'* 

-0.40 

0 



*0 

2.30 X 10'* 

-0.90 

-1700 

12 

OH HO2 ^02 + H2O 


2.90 X 10‘* 

0.00 

-500 

13-^ 

OH -j- H2O2 HO2 -f H2O 


1.75 X 10'^ 

0.00 

320 




5.80 X I0'‘‘ 

0.00 

9560 

14^ 

HO2 4 " HO2 ^ O2 4 * H2O2 


1.30 X 10" 

0.00 

-1630 




4.20 X I0''‘ 

0.00 

12000 

15^ 

0 + 04 *M^ 02 -I-M 


1.20 X 10'’ 

-1.00 

0 

16^^ 

O-FH-FM^OH-FM 


5.00 X 10'’ 

-1.00 

0 

17' 

H2O -F M 


2.20 X \(P 

-2.00 

0 

18-'’ 

H4*H-FM:?i:H2-FM 


1.00 X 10'* 

-1.00 

0 


Note: species are N2* O2, H2, H2O, OH, H, O, HO2, H2O2. 

^AII reactions are reversible. 

^kf{T) = AT*^ t\^{-E^cx/RT)\ units are in cal, mol, cm^, and s. 

does not include 02» H2O, or N2; all collision efficiencies = 1.0. 
^Collision efficiencies for M; H2 = 2.0, H2O = 6.0, all others = 1.0. 
^Troe parameters: a = 0.7346, T*** = 94, T* = 1756, T** =5182. 
^ Rate coefficient is non- Arrhenius: it 13 = it^ + it^ 

^Rate coefficient is non-Arrhenius: k\^ — kc + k^ 

^Collision efficiencies for M; H2 = 2.4, H2O = 15.4, all others = 1.0. 
'Collision efficiencies for M; H2 = 0.73, H2O = 3.65, all others = 1.0. 
•^ Collision efficiencies for M; H2 = 1.7, H2O = 7.0, all others = 1 .0. 


Numerical Method 


The governing quasi 1-D Euler equations are solved 
in finite-volume form throughout the entire domain. The 
equations are solved in split form, in which a fluid convec- 
tion subroutine solves Eq. 1 without the chemistry source 
vector W, 


s ajc 


(13) 


and a finite-rate chemistry integration subroutine solves 
Eq. 1 as a system of ordinary differential equations ignor- 
ing F and H, 



(14) 


The code utilizes the Strang® second-order time step 
splitting approach to couple the fluid convection and finite- 
rate chemistry subroutines. For each complete time step, 


the chemistry subroutine is first called for one-half of^a 
time step, followed by a full fluid convection time step, and 
then followed by another half-step of the chemistry routine. 
Thus, in operator notation, the solution U at time-step n H- 1 
is 

U"+' = (15) 

where £/ represents the fluid subroutine solving Eq. 13 
and Cc represents the finite-rate chemistry subroutine solv- 
ing Eq. 14. The time step splitting approach allows both 
high-quality fluid and chemical solvers to be developed and 
tested independently, and then joined together in relatively 
straightforward fashion. As described by Oran and Boris, ^ 
time-step splitting works well when relatively small time- 
steps are used. Since the global time-step in the present 
code is strictly limited by the explicit Courant-Friedrichs- 
Lewy (CFL) condition, this is not thought to be a problem. 
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The fluid solver used here is the explicit, 2nd-order ac- 
curate (in time and space), symmetric-TVD algorithm de- 
scribed by Yee.*^ The solver employs Roe’s approximate 
Riemann solver modified for nonequilibrium ideal gases. ^ * 
It also incorporates suggestions by Larrouturou^^ to ensure 
species positivity. 

As with most reactive flow problems, the time- 
integration of the chemistry mechanism in this work re- 
quires a stiff-ODE solver to ensure accuracy. The method 
used here employs a semi-implicit trapezoidal scheme as a 
predictor. 





(U"+' - U”) = AtW" 


(16) 


where I is the identity matrix, J is the Jacobian of W, J = 
9W/9U, and is the predicted next value of U. Newton 
iteration of the trapezoidal scheme is used as a corrector, 




(U"+' - U"+') = 


[U" - U"+‘ + ^Af(W" + W"+')] 


(17) 


where is the corrected next value of U. 

The chemistry subroutine is designed to use several lo- 
cal time steps, if necessary for accuracy, within a given 
global half-step. The relative magnitude of the corrector 
(U«+i _ is used as basis for adjusting the 

local time step, similar to the accuracy-based time step 
procedure described by Mott.*^ The Jacobian J is evalu- 
ated analytically once per local time step at U”. Thus, the 
LU decomposition of the bracketed matrix terms on the left 
hand side of Eq. 16 can be re-used in Eq. 17. 


PDRE Simulations 

The PDRE system studied here (Fig.2) is highly ideal- 
ized, consisting of a constant- area (2.0 cm diameter) det- 
onation tube 16cm in length. One end is closed and the 
other end open to the environment, or attached to a simple 
conical converging-diverging nozzle section. The nozzle 
initially converges at a constant 14® angle to a throat 1 .8 cm 
in diameter, and then diverges again at a constant 14® to a 
final exit diameter specified by the user. The computational 
domain automatically scales to fit the nozzle geometry. The 
detonation tube is pre-filled with a gaseous propellant mix- 
ture with no initial velocity. Stoichiometric H 2 -O 2 at an 
initial pressure of 1 atm, and initial temperature of 300 K 
is utilized for all calculations in this work. An idealized, 
massless diaphragm isolates the propellant mixture from 
the nozzle and ambient environment until ruptured by the 
detonation wave. The nozzle section is initially filled with 
H 2 gas at the specified ambient pressure. 

The detonation in each simulation is initiated by specify- 
ing equilibrium combustion conditions (computed at con- 
stant p and €) in the cell adjacent to the closed end-wall. 

Ghost cells are utilized to specify the boundary condi- 
tions in the problem. A reflection-type boundary condition 


Idealized 

diaphragm 


2 



Converging-diverging nozzle 

Fig. 2 Schematic of PDRE geometry studied. 


is utilized at the closed end of the detonation tube to sim- 
ulate a solid wall. The method of characteristics*^ is used 
to calculate the exit flow boundary condition. For sonic 
(choked) or supersonic exit flow, all exit flow properties 
are determined by the interior flow. For subsonic exit flow, 
the ambient pressure is specified and the method of charac- 
teristics is used to compute the remaining flow properties. 
If reverse flow is detected at the exit, a subsonic inflow 
boundary condition is calculated using the ambient H 2 gas 
as a constant pressure and enthalpy reservoir. 

Two additional special restrictions are imposed in the 
simulations. In order to simulate the effect of an ideal- 
ized diaphragm, only the detonation tube portion of the 
domain (from the closed wall to the diaphragm location) 
is computed initially. A reflection-type boundary condition 
is specified at the diaphragm location until the pressure in 
the adjacent cell rises 1 .0% above the initial fill value. This 
special restriction is subsequently removed, and the entire 
domain is computed. Additionally, there is a check per- 
formed when the exit flow is supersonic. Since in this case 
the exit boundary conditions are entirely calculated from 
the interior flow, there is no way for the exit flow to return 
to a subsonic condition.*^ Therefore, at each time step a 
check is made to determine if the pressure from standing 
normal shock at the exit is less than the ambient pressure. 
If true, then the normal shock properties are specified in 
the last interior cell, and a subsonic outflow boundary con- 
dition computed at the exit. 

Grid Resolution 

Grid resolution is a critical issue in all finite-rate chem- 
istry CFD models. Due to the extremely fast kinetics of the 
H 2 -O 2 system, numerical modeling of the reaction zone 
of undiluted H 2 -O 2 Chapman-Jouguet (C-J) detonations 
requires an exceptionally fine grid spacing. The C-J det- 
onation velocity of stoichiometric H 2 -O 2 at and initial con- 
dition of 1 atm and 300 K is calculated by the NASA CEA 
thermochemical code*^ to be 2837.1 m/s. Using the calcu- 
lated frozen shock properties at this velocity, and the chem- 
ical kinetics mechanism described previously, the length of 
the C-J reaction zone for this mixture can be estimated as 
~ 0.02 mm. Thus, a grid spacing of roughly 1 /xm would 
be required to resolve the reaction zone. This is not cur- 
rently practical, nor is it necessary for the performance 
estimates of this study. The chief requirement here is that 
the grid resolution accurately calculate the C-J state and 
velocity. 

A grid resolution comparison of pressure and tempera- 
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Fig. 3 Grid resolution comparison showing pressure (upper 
panel) and temperature (lower panel) profiles of a propagat- 
ing C-J detonation wave at f = 30 fis after initiation. 


ture profiles of a propagating C-J detonation, t — 30 fxs 
after initiation, is shown in Fig. 3. The resolutions vary 
over one order of magnitude. As even the finest resolution 
here, Ax = 0.04 mm, is twice the estimated reaction zone 
size, none of the resolutions can accurately capture the ig- 
nition process behind the leading shock wave. Rather, in 
all cases, the propagating C-J detonation is smeared over 
several cells, resulting in the pressure spike and subsequent 
relaxation toward the C-J state in the finer grids. All four 
resolutions capture the detonation velocity with reasonable 
accuracy (~ 2836 m/s for Ax = 0.04 mm, ~ 2840 m/s for 
the other resolutions), though the precision of this estimate 
is naturally lower for the coarser grids. Note, however, that 
in two finer grids the pressure relaxes to the C-J value in 
a significantly shorter distance than in the coarser grids. 
Based on these observations, a uniform grid spacing of Ax 
= 0.10 mm was utilized for all simulations in this study. 

Thrust Calculations 

The time-dependent thrust is calculated at each time step 
by two different methods. One measure of the thrust as- 
sumes a control volume tightly bounding the solid surfaces 
of the PDRE, and is determined by integrating the pressure 
difference across all surfaces. A second measure assumes a 
rectangular control volume encapsulating the PDRE. This 


measure of thrust is determined from the sum of the time 
rate of change of the internal momentum integral across the 
domain, the momentum flux from the nozzle section, and 
the pressure difference across the control volume. In gen- 
eral, there is excellent agreement between the two thrust 
calculations, and the time-integrated impulse calculations 
agree to within 0.1%. In all simulations, the calculation 
proceeds until the pressure at the closed end-wall is equal 
to the ambient pressure. Thus, these simulations should be 
thought of as single-shot results. 

Results and Discussion 

The blowdown pressure ratio is one of the most critical 
factors governing the performance of any rocket-type sys- 
tem. In this study, this parameter is defined as the ratio of 
the initial fill pressure in the detonation tube to the ambient 
pressure. The quasi 1-D CFD model is used to determine 
the optimum expansion ratio, £, for a converging-diverging 
nozzle for a range of blowdown pressure ratios ranging 
from 1-1000. At each pressure ratio, a variety of expansion 
ratios were tested, and a manual search was performed to 
determine the optimum expansion ratio. In each case, the 
optimum expansion ratio corresponds to the maximum to- 
tal impulse. Two example cases will be examined in detail, 
followed by a comparison of the results among several sys- 
tems over the entire pressure range. 

PDRE Thrust and Pressure Histories 

The optimum nozzle exit diameter for the model PDRE 
system at a blowdown pressure ratio of 1 is found to be 
2.2 cm. This corresponds to an expansion ratio of 1.49. 
Since the manual search for optimal exit radius was con- 
ducted in 1 mm increments, there is a variance on this 
expansion ratio of ±0.26. The mixture-based specific im- 
pulse of the optimized C-D system is 197.48 s. This is 
a modest improvement over the 192.19 s provided by the 
baseline detonation tube without any nozzle at this pres- 
sure ratio. A comparison between the single-shot thrust 
history for the optimized converging-diverging nozzle and 
the baseline detonation tube is shown in the upper panel of 
Fig. 4. The corresponding exit pressure history is shown in 
the lower panel of the figure. It is evident from the thrust 
history that the reduced diameter of the throat in the C- 
D system increases the overall blowdown time compared 
to the baseline detonation tube. Additionally, the thrust 
for the optimized C-D system is actually lower than that 
of the baseline detonation tube for a significant portion of 
the early blowdown history. The main benefit of the C-D 
nozzle is derived from later in the blowdown history when 
the supersonic exit flow from the nozzle provides slightly 
greater thrust than the choked flow of the baseline tube. 
Note that flow from the C-D nozzle becomes overexpanded 
late in the blowdown history, and eventually a normal shock 
forms at the exit by the procedure described in the previous 
section. There is even a small period of negative thrust at 
the very end of the history, as shown by the peak impulse 
marker in the upper panel of the figure. It is evident that 
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Fig. 4 Comparison of thrust (upper panel) and exit pressure 
(lower panel) histories for PDRE systems at a blowdown pres- 
sure ratio of 1. Propellant mixture: stoichiometric H 2 -O 2 . 
Propellant initial conditions: Pinitial ^ ^ '^initial = ^ 

an optimized fixed C-D nozzle can provide only marginal 
benefit 3% additional specific impulse) over the base- 
line detonation tube at a blowdown pressure ratio of 1. 

In contrast, substantial benefits can be realized from a C- 
D nozzle at larger blowdown pressure ratios. The optimum 
nozzle exit diameter for the model PDRE system at a blow- 
down pressure ratio of 1000 is found to be 18,2 ± 0.2 cm. 
This corresponds to an expansion ratio of 102 ± 2. The 
mixture-based specific impulse of the optimized C-D sys- 
tem is 400. 1 1 s, a significant improvement over the 262.04 s 
provided by the baseline detonation tube at this pressure ra- 
tio. A comparison between the single-shot thrust history for 
the optimized converging-diverging nozzle and the baseline 
detonation tube is shown in the upper panel of Fig. 5. The 
corresponding exit pressure history is shown in the lower 
panel of the figure. As is evident from the figure, the blow- 
down process (to 0.001 atm) is considerably longer in this 
case than in the previous one. Note that the optimized C- 
D nozzle provides considerably higher thrust, compared to 
the baseline tube, throughout the blowdown history. This is 
due to the familiar result from classical compressible flow 
that the best performance from a supersonic rocket nozzle 
is obtained when the pressure at the exit plane is expanded 
to the ambient value. As may be seen in the exit pres- 
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Fig. 5 Comparison of thrust (upper panel) and exit pressure 
(lower panel) histories for PDRE systems at a blowdown pres- 
sure ratio of 1000. Propellant mixture: stoichiometric H 2 -O 2 . 
Propellant initial conditions: Pinitial = I Tjniija] = 300K. 

sure history, the C-D nozzle expands the exhaust flow by 
roughly three orders of magnitude compared to the base- 
line tube. Thus, while the fixed nozzle cannot dynamically 
adapt to provide optimum expansion throughout the entire 
blowdown process, the resultant performance gain is still 
quite significant (^ 53%). It is worthwhile to note that the 
peak impulse (as shown by the marker in the upper panel) is 
reached by 1 141/xs, roughly half the entire blowdown time. 
After this point, the thrust oscillates both slightly above and 
below zero, leading to a small loss in impulse by the end of 
the blowdown process ('^ 0.4%). This suggests that, in an 
actual cyclic PDRE system using this geometry, purging 
and propellant refill should begin much earlier than allow- 
ing the end- wall pressure to decay to the ambient pressure 
value. 

Equivalent SSRE Model 

It is instructive to compare the specific impulse of both 
PDRE systems with a SSRE under equivalent modeling 
assumptions. While frozen and equilibrium rocket perfor- 
mance calculations can be readily obtained from the CEA 
code, these results are not directly comparable to the finite- 
rate chemistry model used in the current PDRE code. The 
primary concern is the tendency for chemistry to slow down 
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in real nozzle systems as the temperature and pressure are 
reduced in the expansion process. Thus, it is best to com- 
pare the finite-rate PDRE systems with a finite-rate SSRE 
model. This SSRE CFD model is heavily derived from the 
PDRE code. The 16 cm detonation tube is replaced with 
a thrust chamber 0.1 cm in length, and 3.6 cm in diameter 
(Fig. 6). This thrust chamber then converges at a constant 
14® angle to a throat 1.8 cm in diameter. After the throat, 
the nozzle again expands at 14® to the exit diameter speci- 
fied by the user. The throat diameter and nozzle expansion 
rate are thus identical in both the PDRE and SSRE models. 
The equilibrium temperature, pressure and composition of 
stoichiometric H 2 - 02 , burned at constant pressure and en- 
thalpy from initial conditions of 300 K and 1 atm, are fed 
as a constant enthalpy reservoir inflow boundary condition 
into the domain. 

Similar to the PDRE nozzle optimization study, the 
finite-rate SSRE CFD model is run at various expansion 
ratios until the optimum specific impulse is obtained for a 
given pressure ratio. However, as the solution converges 
to steady-state in this model, flowfield information from 
the highest pressure ratio case can be used to guide area 
selection at the lower pressure ratios. In each case, the 
SSRE model is run until the solution converges. In general, 
the specific impulse and optimum expansion ratio using 
finite-rate chemistry are slightly larger than that for frozen 
chemistry. Additionally, if the chemistry is frozen in the 
SSRE CFD model, there is excellent agreement (to within 
0.3% in specific impulse) with the frozen-chemistry predic- 
tions of CEA. 

Performance Comparison 

The mixture-based specific impulse for the baseline det- 
onation tube, the optimized fixed C-D nozzle PDRE sys- 
tem, and the equivalent optimized SSRE system is plotted 
for blowdown pressure ratios of I, 2, 5, 10, 20, 50, 100, 
200, 500 and 1000 in the upper panel of Fig. 7. The cor- 
responding area ratios for the C-D systems are shown in 
the lower panel of the figure. As would be expected from 
the example results discussed previously, the relative gain 
from the C-D nozzle systems becomes more pronounced 
at higher pressure ratios. It is interesting to note that that 
both the baseline and C-D nozzle PDRE systems outper- 
form a SSRE at pressure ratios below ~ 7. This is due 
to the fact that the relative gain from the extra pressuriza- 
tion provided by the detonation process is greatest in this 
regime. Thus, the PDRE may have considerable theoreti- 




Fig. 7 Performance (upper panel) and area ratio (lower 
panel) comparison of various pulse detonation and steady- 
state rocket devices. All results obtained using finite-rate 
quasi 1-D CFD calculations. The final blowdown pressure at 
the end-wall in the PDRE is equal to the ambient pressure 
(Pflnal = Pambient)* Propellant mixture: stoichiometric H 2 -O 2 . 
Propellant initial conditions: Pinitial = ^ '^'initial = 

cal potential for rocket-type applications when the pressure 
of the ambient environment is high. Additionally, at higher 
blowdown pressure ratios, a PDRE with an optimized, fixed 
C-D nozzle has a greater specific impulse than an equiva- 
lent SSRE. These results indicates that a PDRE equipped 
with a fixed C-D nozzle, though obviously not as efficient 
at generating thrust as a (theoretical) dynamically adaptive 
nozzle, can still yield performance superior to an equivalent 
SSRE over a wide range of pressure ratios. Note, how- 
ever, that this performance gain becomes relatively smaller 
at higher pressure ratios. 

Summary and Future Work 

A quasi 1-D, finite-rate chemistry CFD model for study- 
ing PDRE gasdynamics and performance is described and 
implemented. The performance of a simple fixed, but op- 
timized, converging-diverging nozzle design is compared 
with a baseline detonation tube over a range of blowdown 
pressure ratios from I-IOOG. The results demonstrate that 
even relatively simple fixed nozzle designs can make sig- 
nificant improvements in PDRE performance at high pres- 
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sure ratios. Additionally, as practical systems would likely 
employ cryogenically-stored, fuel-rich propellant mixtures, 
work is planned to repeat the optimization study for a range 
of stoichiometries and initial temperatures. Work is also 
underway to implement an ethylene chemistry model for 
comparison with experimental results from the literature. 
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